Adiabatic pumping through interacting quantum dots 
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We present a general formalism to study adiabatic pumping through interacting quantum dots. 
We derive a formula that relates the pumped charge to the local, instantaneous Green's function of 
the dot. This formula is then applied to the infinite-[/ Anderson model both for weak and strong 
tunnel-coupling strengths. 
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Introduction. The idea of producing a DC current at 
zero bias voltage by changing some parameters of a con- 
ductor periodically in time dates back to the work of 
Thouless . This method of exploiting the explicit time 
dependence of the Hamiltonian of the system is known 
as pumping. If the parameters change slowly as com- 
pared to all internal time scales of the system, pumping 
is adiabatic, and the average transmitted charge does not 
depend on the detailed time dependence of the param- 
eters. For non-interacting mesoscopic systems, Brouwer 
@ , using the concept of emissivity proposed by Biittiker 
et al. 3[, related the charge pumped in a period to the 
derivatives of the instantaneous scattering matrix of the 
conductor with respect to the time-varying parameters. 
In case of noninteracting electrons, a general framework 
for the computation of the pumped charge has been de- 
veloped 0, 3 • Pumping through open quantum dots has 
also been investigated experimentally 

The situation is profoundly different for pumping 
through interacting systems. In fact.there are only few 
works that address this problem |E 13 with methods 
suited to tackle specific systems or regimes. As far as 
pumping through interacting quantum dots is concerned, 
the work by Aono Q exploits the zero-temperature map- 
ping of the Kondo problem to a noninteracting system 
and uses the noninteracting formalism. On the other 
hand, Cota et al. 8] study adiabatic pumping in a 
double-dot system in the sequential tunneling limit. 

The aim of this Letter is to derive a formula for the 
charge pumped through an interacting quantum dot, 
which is valid from the high-temperature limit where se- 
quential tunneling dominates down to low temperatures 
where Kondo correlations are relevant. 

Model and formalism. We consider a single-level quan- 
tum dot coupled to two noninteracting leads. The system 
is described by the Hamiltonian 



H = H\ 



leads 



H. 



dot 



H t , 



(1) 



With i/lcads = J2k,a.a 6 a( k ) c lka C ^a, where C ak a (c* fca ) 

is the fermionic annihilation (creation) operator for an 
electron with spin a =1,1 and momentum k in lead a — 
L, R. The leads are assumed to be in thermal equilibrium 



with the same chemical potential and to have flat bands 
with constant density of states p a . 

The quantum dot is described by iJdot = [e + 
Ae(i)] Yla n cr + Un-^ny with n a = d) a d a , where d a (dj.) is 
the fermionic annihilation (creation) operator for a dot 
electron with spin a. The level position of the dot con- 
tains a time-independent part e and a time-dependent 
part, Ae(t). Coulomb interaction in the dot is described 
by the on-site energy U. Tunneling is modeled by H tun — 



H.c. 



with time-dependent tunnel 



matrix elements V a {t). We only allow for the modulus, 
but not the phase, of V a (t) to vary in time, since a time- 
dependent phase would correspond to a bias voltage. 

By periodically changing (at least two of) the three 
quantities VL(i), Vn(t), and Ae(i), a finite charge can be 
pumped through the quantum dot. The charge Q that 
is pumped after one cycle T is connected to the time- 
dependent current Jh(t) flowing through the left barrier 
via the relation Q = J J^(T)dT. The starting point 
for our analysis is the exact relation that expresses the 
current in terms of the dot Green's function [lfj 



Mt) = -t 



2e 

h 



Im 



r L (t,t) 



dt'e- iw{ - t '- t ^ n Y lj {t',t)G t rTrT {t,t') 



(2) 



with r L (ti,i) = 2Trp L V h (t)V{{t 1 ), and f(oj) is the Fermi 
function. The lesser, retarded, and advanced Green's 
function are defined as usual, G< a (t,t r ) = i(d' a (t')d cr (t)), 

Gl a {t,t') = -i0(t-f)<K(t),4(*')}>. and G^(M') = 
[GJ J(T (t', t))* . The Green's functions are diagonal in spin 
space since tunneling is spin conserving. Furthermore, 
spin degeneracy yields G^(t,t') — Gn(t,t') = G(t,t'). 
We remark that the Green's functions G(t, t') are de- 
fined with a Hamiltonian that explicitly depends on time. 
They are determined by the Dyson equation 

G(t,t') = g[t,t')+ [ *idt 2 (5(Mi)S(ti 1 ia)S(*2,t / ). ( 3 ) 



in matrix notation A 



A T A< 
A a 



for the bare, g, and 



full Green's function G, and the self-energy E. The lat- 



2 



ter takes into account the tunnel coupling T a (t) to the 
leads, the Coulomb interaction U in the dot, and the time- 
dependent part Ae(t) of the level position. Note that 
*2) = £(*i>£ 2 > {H( T )} T £[ti t 2 ]> IS a functional of the 
time-dependent Hamiltonian H(t) on the interval \t\, £ 2 ]. 

We are interested in the behavior of the self-energy 
and, thus, the Green's function for a slowly varying 
Hamiltonian H(t). This means that the time scale 
over which the system parameters are varying is large 
compared to the lifetime of the system. To construct 
the adiabatic expansion of the self-energy we first lin- 
earize the time dependence of the Hamiltonian, if (r) — > 
if (to) + ( T — to)H(t ), with respect to some fixed time to, 
and expand the self-energy up to linear order in the time 
derivative, where the time ordering in if (to) is still done 
with respect to time r. The relation Jjj 2 TH{to)dr — 

(ti + r 2 )/2 H(to)dr, valid for each segment of time 
evolution between two vertices at times t\ and r 2 in 
the self-energy, motivates a global replacement of the 
time variable r with the average time (t± + £ 2 )/2 in 
the self-energy. This replacement defines an approxi- 
mation, which we refer to as the average-time approx- 
imation^^. As a result, the dependence of the self- 
energy on the function H(t) over the interval [£1, ^2] is 
replaced by the dependence on the three times to, t%, 
and £2 only, and we arrive at the adiabatic expansion 
S(£i,£ 2 , {i£(r)} re[ti t2] ) -> S (£i,£ 2 ,£o) + Si(£i,£ 2 ,£o) 
with 



G<(w,t) 



-2if(u>)1mG c (w,t), where G (w,£) 



E (£i,£ 2 ,£o) = S(£i,£ 2 ,{if(£ )}) 
Si(£i,£ 2 ,£o) = 



(4) 



£i+£ 2 \ 9E (£i,£ 2 ,£o) 
— Wo ' (5) 

The lowest term in the adiabatic expansion corresponds 
to replacing the time-dependent Hamiltonian H (r) with 
the constant value ff(£o). Then, £o(£i, £ 2 , £0) depends 
on £1 and £ 2 only via the difference £1 — £ 2 , and we 
can introduce the Fourier transform Eo(^,£o) = J d(ti — 
£ 2 ) exp[«o;(£i - i 2 )/fi]£ (ti, £ 2 , £0). 

The adiabatic expansion G(£, £') — > Gq(£, £', £0) + 
Gi(£, £', £0) for the Green's function follows from that for 
the self-energy via the Dyson equation Eq. . Again, we 
can introduce Fourier transforms Go/i(w,£o) = / d(t — 
£') exp[iw(£ — £')/7L]Go/i (£,£', £o)- Since our goal is an 
adiabatic expansion of the current at time £ as given in 
Eq. J5J, we choose from now on £q = £. This results in 



GoM) = -t (uj,t) 



r 1 A . fc aG (^,£) dt (uj,t) ^ 
Gi(w,£) = in 5 ^ G (w,£) 



(6) 



duo 



Of 



ifc*, , .d 2 lbo(u,t)- 

We specify these matrix equations for the retarded 
and lesser part and make use of the equilib- 



[u> — e — S (cj, £)] . Furthermore, the adiabatic expan- 
sion for Fl(£', £) can be constructed as Fl(£', £) — ► Fl(£) — 
^r L (£) with r L (£) = T L (£, t). Plugging everything into 
Eq. J2J we find that the zeroth-order term of the adi- 
abatic expansion for the current vanishes (as it should 
since it is equivalent to time-independent problem at 
equilibrium). The first-order correction is given by 



■h(t) = du[- 



duj 



Rc 



dt 



[r L (£)G ( W ,£)] 



(GK^-'G^t) 



(8) 



A factor 2 accounts for the spin degeneracy. Equation (jHJ) 
is the central result of this Letter ^2]. It generalizes 
Brouwer's formula to interacting quantum dots. We 
emphasize that this result relies on the average-time ap- 
proximation for the self-energy. The latter is exact when- 
ever the self-energy contains two vertices (either tunnel- 
ing or interaction) only. This is the case for U = 
but also for U — > 00 as long as the self-energy is cal- 
culated up to linear order in the tunneling coupling T, 
as well as for arbitrary interaction at zero temperature, 
where the interacting problem can be mapped to a non- 
interacting one. We now specialize Eq. (JHJ) to the case 
of weak pumping due to time-dependent tunneling bar- 
riers, T Q (£) = f a + AT a (t), where |Ar Q (£)| < f a at 
any time £. To the lowest order in Ar Q (£) the charge 
Q = J Ji^{r)dT in one period T is 



Q = - 



er/T 
7rf L f r 



1 "ajj ST rM ' 



(9) 



where T = _T L + Ir, and n = J Q T AT L {t) AT K {t)dt. 
The symbol S(ui) denotes the phase of the Green's func- 
tion Gq(w) = |G (w)| exp[iS(uj)] computed with T a (t) re- 
placed by f Q , and T(oj) = 2f L f R /f ■ Im[G Q (w)] can be 
interpreted as the transmission probability through an 
interacting quantum dot 

Examples. We now consider only weak pumping 
with the barriers, and we restrict ourselves to the case 
Fl = Tr = r/2. We start by studying the nonin- 
teracting single- level quantum dot, using Eq. J5J, with 
Gq(w) = (lu - e + |r From inspection of Eq. © it is 
clear that there is no pumping in the noninteracting case 
if the level is resonant (e = 0). In the high-temperature 
limit j3T <C 1 (with (3 = I/ZcrT 1 ), the pumped charge 
reads Q = -f /"(e). 

We now turn our attention to weak pumping with the 
barriers in the limit of large electron-electron interaction, 
U — > 00. For temperatures larger than the Kondo tem- 
perature (defined below), we approximate the instanta- 
neous Green's function of the dot within the equation-of- 
motion method [FJ|. Replacing r a (£) by f Q one finds 

G H = (i-(«))(^^^H + # + /H) _I , 



rium relations T^Q-(u),t) 



-2i/(w)Im£&(w,t) and 



-10 



o_ 

e/r 
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FIG. 1: Pumped charge in units of erj/T 2 as a function of the 
level position in units of f for U — > oo (solid line) and U — 
(dashed line). The temperature is Ub.T = 5f , and E c = 20f . 



where = ± {v (§ + ^) - Rc 

and (n) = — J ^rlm |Go(w)} f(u>) is the occupation of 
the level per spin, ip the Digamma function, and E c 
a high-energy cutoff. For the high-temperature limit, 
(3T <C 1, we obtain the analytical expression 



erj 
T 



/"(e) + 



/'(e) 



2A(e)/f 



(10) 



/(e) V" w l + /(e) 

Figure ^ shows the pumped charge as a function of 
the level position in the high-temperature limit. The 
enhancement of the pumped charge as compared to the 
nonintcracting case, is mainly due to the fact that, in 
the presence of interactions, the bare level is renormal- 
ized by an amount which depends on T(t), and hence 
the level position becomes time dependent. The oscil- 
lation of the level increases the pump effect [third term 
in Eq. {TP))] (lif . Also the fact that the level width is 
energy dependent, E[l + f(u>)], has some small effect on 
the pumped charge [second term in Eq. (fTUf) ] . We note 
that for j3 — > 0, the third term in Eq. O goes as (3 2 f E c , 
while the other two terms go as (/3f ) 2 16]. The shape of 
the curves in Fig. ^ are easily understood from the de- 
pendence of 85 /8T around \oj\ < k B T on the bare level 
position e. In the nonintcracting case, the scale on which 
the phase 5 varies around the level position e increases 
linearly with f ; i.e., dS/dT changes sign when tuning the 
level position through the Fermi energy. In presence of 
interaction, though, the dominant mechanism is the vari- 
ation of the level renormalization, which shifts 5(lu) along 
the uj axis, with no sign change in OS/dT. 

The equation-of-motion method gives qualitative, re- 
liable results down to the Kondo temperature, given by 
k B T « k B T K = V r ^T/2exp(-7r|e|/f). In Fig. EJ we 
show the temperature dependence of the pumped charge 
for the interacting quantum dot, obtained by numerical 
integration of Eq. @ and for comparison, the noninter- 
acting result. At very high temperatures the pumped 
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(k B T)/|8| 

FIG. 2: Pumped charge in units of eq/T 2 as a function of 
temperature (in units of the level position), for U — ► oo (solid 
line) and for (7 = (dashed line). The level position is fixed 
at e = — 2.5F. In both cases a maximum appears at almost 
the same temperature. The cutoff energy is E c = 20r. 
Inset: Pumped charge as a function of temperature (in units 
of Ik) for U -> oo. For < T < 0.5T K , Q is obtained by 
means of the mean-field slave-boson method , performing a 
Sommerfeld expansion. For T > 1.5Tk it is computed numer- 
ically using the equation-of-motion Green's function. 



charge tends to zero. Decreasing the temperature the 
charge exhibits a maximum both in the interacting and 
nonintcracting case. It occurs when the level position 
and the temperature are of the same order. Its position 
is determined by the spectral weight of the integrand 
function in Eq. J^J which falls in the energy window set 
by temperature through the derivative of the Fermi func- 
tion. Approaching the Kondo temperature, the pumped 
charge in the interacting system increases rapidly, indi- 
cating Kondo correlations. 

To address the limit T <C Tk we resort to the slave- 
boson method [TtI ] in the mean-field approximation in 

the boson field. The instantaneous dot Green's func- 

r ft) 

tion can be written as G T (ui,t) = r(t) ^pf^' ^' wnere 
the pseudofermion Green's function is Gp f (u;,t) = (u> — 
e p f(t) + ir p f(t)/2) _1 . We are interested in G pf (w), where 
r Q (t) is replaced by r Q . The renormalized level position 
e p f and the renormalized rate f p f have to be found as the 
solutions of the non-linear system of equations: 



du> 



2f / _Re{G pf H}/H + 6 



2tt 
2^' 



pf 



Im{G pf M}/H = l 



e = 0, (11a) 
(lib) 



IV 
r 



At zero temperature the pumped charge can be expressed 
by means of Friedel's sum rule and Eq. (JSJ, as 



Aer/ d(n) 



■ sin 



(n(n)) , 



(12) 



r or 

which relates the pumped charge to the average occupa- 
tion per spin (n) only. The full knowledge of the latter, 
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FIG. 3: Pumped charge in units of en/T 2 as a function of the 
level position in units of F, for U — > oo, T — and E c = 20F, 
computed by means of the mean-field slave-boson method. 

e.g. from numerical renormalization group, would estab- 
lish an exact solution of the problem. Within the mean- 
field slave-boson approach we get (n) — 1/2(1 — f p f/f) 
fl9|. and F p f is computed from Eqs. (|TT|) . In the unitary 
limit (-e > f pf and T < T K , such that (n) -> 1/2) 
the level is renormalized to resonance e p f — > 0, and the 
pumped charge is zero. This result is consistent with the 
fact that in the unitary limit the problem maps to the 
noninteracting dot with the level shifted to resonance, 
and that for the free-electron case there is no pumped 
charge when the level is at the Fermi energy. On the other 
hand, in experimentally relevant situations the renormal- 
ized level is not exactly at the Fermi energy and non- 
negligible charge pumping occurs. 

In Fig. we show the charge pumped at zero tempera- 
ture obtained solving numerically Eqs. (fill) . As expected 
the charge tends to zero when the level is deep enough 
below the Fermi energy. The behavior of the pumped 
charge around T <C Tk can be obtained by performing a 
Sommerfeld expansion in Eq. @, and in Eqs. i|ll[l . The 
pumped charge goes as T 2 . Comparing the temperature 
behavior of the charge for T <C Tk with the one for T 
just above Tk [see inset of Fig. EJ, we expect a maxi- 
mum at around Tk- Roughly speaking, this extremum is 
analogous to the one that occurs at higher temperatures. 
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